A dynamical evolution model on the black hole horizon 
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This paper demonstrates a dynamical evolution model of the black hole (BH) horizon. The result 
indicates that a kinetic area-cells model of the BH's horizon can model the evolution of BH due to 
the Hawking radiation, and this area-cell system can be considered as an interacting geometrical 
particle system. Thus the evolution turns into a problem of statistical physics. In the present work, 
this problem is treated in the framework of non-equilibrium statistics. It is proposed that each area- 
cell possesses the energy like a microscopic black hole, and has the gravitational interaction with the 
other area-cells. We consider both a non-interaction ideal system, and a system with small nearest- 
neighbor interactions, and obtain an analytic expression of the expected value of the horizon area 
' of a dynamical BH. We find that, after a long enough evolution, a dynamical BH with the Hawking 

radiation can be in equilibrium with a finite temperature radiation field. However, we also find that, 
the system has a critical point, and when the temperature of the radiation field surrounding the BH 
approaches the critical temperature of the BH, a critical slowing down phenomenon occurs. 

j^j . PACS numbers: 04.70.Dy, 04.60.Pp 

The discovery of the thermal radiation of black holes (the Hawking radiation [1] ) is a significant breakthrough in 
£T) ■ our understanding of black holes. However, after three decades, there is still disagreement on a number of important 
7-H ' issues. For example, it is commonly accepted that black holes have finite temperature, which is proportional to the 
surface gravity of the event horizon and inversely proportional to the mass, and Black holes can be in equilibrium with 
a finite-temperature radiation field. However, it is also known that the thermal radiation of black holes is a quantum 
effect. As the thermal radiation goes on, the black hole is emitting particles. As its mass decreases, the temperature 
rises and, as a result, the black hole emits particles more quickly. It is thus foreseeable that the black hole cannot 
keep thermal equilibrium, and it will finally evaporate and reach zero mass. Now we turn to review the formation of 
a black hole from the classical point of view. A star in its late stage of evolution may evolve into a black hole, due 
I ' to the gravitational collapse. This black hole quickly reaches a stationary state, which can be characterized by three 
parameters: mass, charge and angular momentum. Even if we consider that the black holes with nonzero angular 
momentum and charge are in an active excited state, the Schwarzschild black holes with zero angular momentum and 
zero charge are in an inactive ground state, which is the end-result of star evolution. In light of the mentioned-above 
O" 1 points, questions naturally arise: Certainly a black hole with a certain temperature can be taken as a thermodynamic 
| \ system, but can it be taken as a system in thermal equilibrium? Is equilibrium statistics applicable? Will thermal 
bJ[). radiation really lead to the total evaporation of a black hole? What is the end-result of star evolution? They are still 
open questions. Before the real theory of quantum gravity is completely established, can we find consistent answers? 

The present work, we aim at the establishment of a model of non-equilibrium dynamical evolution, in the framework 
of statistical physics. The expected area of the event horizon is calculated, and its long time evolution is observed. 
These results are used to shed light on the mentioned-above questions. 

The model. In the past years, Loop Quantum Gravity (LQG) has become a serious candidate for a non-perturbative 
quantum theory of gravity [2]. Its most notable prediction is the quantization of geometry [3]. As an analogy, in 
LQG, the fabric of space is like a weave of tiny threads, and each thread poking through a surface gives it a little bit 
of area. The surface area of a black hole (BH), then, is generated by all the threads puncturing it. The event horizon 
is flat except at these punctures, where it can flex, and the microstates of the BH are defined by the different ways 
the event horizon can flex in or out [4] . 

According to the results of LQG, the horizon surface is split into a series of N discrete cells (small surface elements 
at the Planck scale). Each cell has a discrete spectrum of area 
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Ai(ji) = SnllWjiUi + l), i = l,2, 3, N; j t = 0, 1/2, 1, 3/2, • • • , (1) 

where l p is the Planck length, 7 is the Immirzi parameter, and j is the irreducible representation of the gauge group 
SU (2). Due to the Hawking quantum evaporation, the geometrical area of these microscopic cells evolves with time. 
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Now, we can investigate the dynamical evolution of the event horizon of a BH from the quantum point of view. 
During the time 1 interval of t ~ t + dt, the event horizon of a macroscopic BH can be considered as a set of area- 
cells at the Planck scale. The area-cells, with half-integer spins, j = 0,1/2, 1, 3/2, ••• (irreducible representation 
of the SU(2) gauge group), are characterized by quantized geometrical area, and thus are referred to as quantum 
geometrical particles 2 . For an individual area-cell, it is a small quantum BH, and its energy has the same form as 
a Microscopic Black Hole. However, adjacent area-cells have gravitational interaction, and many area-cells together 
form an interacting system of geometrical particles with spin. Since there can be various possible transitions among 
the area-cells in different states, the collective transition of the system becomes the quantum evaporation of the 
macroscopic BH. We use a heat bath at temperature T to simulate the radiation field surrounding the BH, and 
suppose that the geometrical particle system is in contact with this heat bath, but it may be in a non-equilibrium 
state. 

A system consisting of these geometrical particles with area (or spin) as a characteristic variable is very similar 
to a spin-lattice system in condensed matter physics. The contact of the system with a heat bath may reflect the 
situation of a BH in a radiation field. The Hawking evaporation leads the characteristic parameters to evolve with 
time. In this case, the evolution of the event horizon of a BH can be investigated by studying the dynamical behavior 
of a geometrical particle system in non-equilibrium states. A generalized Glauber-type dynamics with single spin 
transition mechanism [6-8] can be very useful here. In this paper, we will reformulate the Glauber dynamics to make 
it suitable for the kinetic geometrical system, and then apply it to a Black Hole with discrete area spectrum. 

The formulation of the dynamical scheme. The area of each fixed cells on the event horizon of a BH is represented as 
a stochastic function of time Ai(t), (i = 1, • • •, N), which can be taken as discrete values. Transitions can occur among 
these values. The transition probability Wi(Ai(t) — > Ai(t)) from configuration (A\(t), ^(i), • • • , Ai(t), • • • , Ajv(i)) to 
configuration (^A\(t) , A%(t) , • • • , Ai(t), • • • , Aw(tfj , in general, depends on the momentary values of the neighboring 

cells as well as on the influence of the heat bath. For this reason statistical correlations exist between different 
geometrical particles. Therefore, it is necessary to treat the entire TV-particle system as a whole. The evolution of the 
particles' area functions, which describes the evolution of the system, forms a Markov process of N discrete random 
variables with a continuous time variable as argument. 

We introduce a probability distribution function P({A}; t), which denotes the probability of the geometrical particle 
system being in the state of {A} = (A\, ■ ■ ■ , Ai, ■ ■ ■ , An) at time t. Let Wi(Ai — > A{) be the transition probability per 
unit time that the ith particle transits from area Ai to another possible area Aj, while the others remain unchanged. 
Then, on the supposition of single-particle transition, we may write the time derivative of the function P({A};t) as 



d_ 
~dt 



P({A}; t) = J2Y,9 ( A i) [~ W i( A i A i) P (i A h t) + Wi(Ai - Ai)P(- ■ ■ , A u ■ ■ ■ ; t)] , (2) 



A, 



where g (Aij is the degeneracy of state Af. This is a probability equation, in which the first term in the right-hand 
side denotes the decrease of the probability distribution function P({A},t) per unit time due to the transition of 
the particle state from the initial value Ai (i — 1, 2, • • • , N) to various possible final values A i7 and the second term 
denotes the contrary situation. We shall refer to Eq. (2) as the master equation since its solution would contain the 
most complete description available of the system. 

It is the most crucial step, obviously, to determine the transition probabilities, Wi{A,i — > Ai), i = 1, 2, 3, • • • , N, 
before the master equation can be possibly solved. Then, how to determine the transition probabilities? We have 
both mathematical and physical considerations. On the one hand, mathematically, the transition probability must 
be positive definite and normalized; and physically, a thermodynamic system in a slowly varying process must have 
ergodicity and satisfy the detailed balance condition. On the other hand, the transition probabilities of the individual 
particles depend mainly on the momentary values of the neighboring particles as well as on the influence of the heat 



1 With respect to "time" , we must point it out that, although it is true that there exists a notion of time in the 3+1 formulation, 
it is certainly untrue that the notion of time is unambiguous. However, in order to carry on the discussion of the dynamic 
evolution of black hole horizon using non-equilibrium statistical mechanics, here we still give a clear definition of "time" as an 
assumption. The extent to which this assumption is acceptable can be judged from the final result of the research. 

2 As far as I know, the idea of "quantum geometrical particle" was first given by Major and Setter in Ref. [5]. What I especially 
would like to emphasize here is that, except the quantized geometric area (labeled by spin j), the other microscopic details of 
the geometric particles have been averaged out, in a way technically similar to the mean-field approximation used to treat the 
ferromagnetic systems in condensed matter physics. This is an issue worthy of further discussion. 
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bath. Thus, the transition probability from Ai to Ai must depend on the heat Boltzmann factor of the neighboring 
particles. Based on these consideration, the following form of the transition probability 



1 



Wi(Ai -> Ai) = — exp 



-pHefflAi^AH 



(3) 



is a natural choice, where Qi is the normalization factor determined by the normalized condition J2a - 9 (^*) ^»(^» — * 

Ai) = 1 and the summation is taken over the neighboring particles of i. 

Usually, we are only interested in the expected value of a certain function of the area 

(f(A k ))=Y,9(A 1 )---g(A N )f(A k )P({A},t). (4) 

{^4} 

According to the definitions (4) and the master equation (2), the time-evolution equation of (/ (A k )) can be derived 
( see Appendix) 



Jt (f(A k)) 



(f (A k )) +'£g(A 1 )...g (A N ) I J2 9 (Ak) f (i fe ) W k (A k A k ) P({A}, t). 



(5) 



Application to the ideal geometrical particle system. As a first step, we consider an ideal geometrical particle system. 
In this case, the total area of horizon, A(t) = A k (t), can be taken as the parameter. Therefore the time-evolution 
equation (5) can be simplified to the following form, 



d 
~dt 



(A(t)) = (A(t)) + ]T g (A) 9 (A) AW {A - A) j P(A; t). 



(6) 



Because 



]T. g (i) AW(A^A) = 



Ea9(a) 



A exp 



-pHeffU 



Eiff(i)exp [-men (a) 



= (A) C 



and 



J2g(A)P(A;t)=TrP(A;t) = l, 



(7) 



(8) 



we can obtain 



(A(t)) = (A) eq + (A(0))-(A) eq < 



(9) 



The expression (9) means that the event horizon as a whole evolves in the form of exponential decrease with time. 
Obviously, the result of a long-time evolution is (A(oo)) = (A) . From Eq. (7), as long as the discrete spectrum 
expression of A is given, the expected value of the horizon's area, (A) , can be calculated. 

Application to interacting geometrical particle system. In our consideration, an individual area-cell (a quantum 
geometrical particle) is regarded as a microscopic quantum BH, and between neighboring geometrical particles there 
is gravitational interaction. We can adopt as a hypothesis that the energy (mass)-area relation is a power law, and as a 
specific choice, we suppose that its energy (mass) is proportional to the square root of the area (Schwarzschild-type) 3 , 



Ei = m l cx y/ Ai{ni), 



(10) 



3 What is chosen here is the same as Ref. [5] (see formula (5) of this paper). 



3 



where Aj(rij) = Aizl^p/y/ rii (rii + 2) = aa^Jni (rii + 2), ao = AttI^, rii — 2ji = 0, 1, 2, 3, • • •. The gravitational interac- 
tion between neighboring geometrical particle, however, is yet an unknown problem. But, our understanding is that: 
(1) the Planck scale l p is the minimum length in the microscopic structure models of the quantum gravity theory; (2) 
the event horizon is flat except at those punctures generated by all the spin network's edges puncturing it [4]. So, we 
might as well adopt simply a classical gravitational potential as follows 



u(r tf ) 



CO, 



r i,j = °< 



— j^y/Ai(rii)y/Aj(nj), r% t j = NN distance, 



0, otherwise 
Thus, the effective Hamiltonian (with only nearest-neighbouring (NN) interaction) can be written as 
Hqt. = ^e t + ^u(r 13 ) = m^[n l (n l + 2)} 1/4 -u^[n l {n i + 2)n j (rij +2)} 1/4 , 

i i<j i 



(11) 



(12) 



where m and u are the scaling factors, and m > 0, u > 0, and each sum runs all NN pair. Thus, the transition 
probability per unit time, that the characteristic quantity of the ith particle transits from one value rii to another 
possible value hi, rii — > hi, can be written as 

Wi (Ai (m) -> Ai (nj)) = ^ cxp j-pro [h t (hi + 2)] 1/4 + [3uJ2 [«< («i + 2) (n J+M , + 2)] 1/4 | , (13) 

where 

Q 4 = ^ ("« + ex P ex P ^ ~/5 m ["» ("» + 2 )] 1/4 + ("« + 2 ) ni+lu ( n *+™ + 2 )] 1/4 f • 

rii I w ) 

Here, the degeneracy of spin state ji, g (ji) = 2% + 1 = hi + 1, is taken into account. At high temperature, the 
summation for n, can be approximately replaced by an integral — ► J dhi. 

In order to get an analytical result, we consider the case of a weak gravitational interaction and a high environmental 
temperature. Then, one can obtain 



^2 9 («») (fit) Wi(n 4 -> fti) = -p-^ao 



+ 2— (n l+w ) + ((—)) 



(14) 



Vj(«,) jA\{h~)Wi(rii -» n,) = 
' v pro 

Therefore, we have the following evolution equations 



1 + --?= S v / ^(«i+») + O ( (-) 
m -,/ao \\mJ 



p2 m 2 m 



equivalcntly, 



20a 

/3 2 m 2 



V G4T(t)) + (v / ^))+-- 
lOv 7 ^ 



16 u 



TV- 1 - 
pm \ pm m 



V^)) + ( v ^)> + . 



(15) 

(16) 
(17) 



(18) 



(19) 
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The exact solution is 



t^eH+I^'^' (20) 

\ /3m m 



where (A(t)) = J2i (Ai(t)) denotes the total area, C\ and Ci are the integral constants, and r is the relaxation time 
of the system, 

We notice that the system has a critical point, f3 c — — — . When the gravitational interaction is negligible, the 
system will evolve rapidly to the equilibrium state, {A)° eq = N j^j^ oc ^iV(^p) 2 , which agrees with the result of 
Ref. [5]. However, a more complex dynamical behavior will appear when the interaction cannot be ignored. If the 
temperature of the radiation field is much lower than the critical temperature T c , the system will still rapidly approach 

the equilibrium state (A) eq — N ^ ?« u - — lj ■ However, if the temperature approaches the critical point, the 

system can hardly reach the equilibrium state, and this phenomenon is commonly known as the critical slowing down. 
It is surprising that the horizon area will increase continuously when the temperature is higher than the critical point. 
But this case should be excluded in the present discussion using this specific method. The reason is that, due to the 
very intensive heat exchange, we cannot regard the radiation field surrounding the BH as a heat bath with constant 
temperature. 

This paper demonstrates a dynamical evolution model of the black hole (BH) horizon with discrete area spectrum. 
The result indicates that the evolution of BH due to the Hawking radiation can be modeled by a kinetic area-cell model 
of the BH's horizon, and this area-cell system can be considered as an interacting geometrical particle system with 
spin. Thus the evolution turns into a problem of statistical physics. In the present work, this problem is treated in the 
framework of non-equilibrium statistics, and the expected area of the event horizon is obtained. We find that, after a 
long enough evolution, a dynamical BH with the Hawking radiation can be in equilibrium with a finite temperature 
radiation field. However, we also find that, the system has a critical point, and when the temperature of the radiation 
field surrounding the BH approaches the critical temperature of the BH, a critical slowing down phenomenon occurs. 

Of course, the present work is only a preliminary attempt on the evolution of the BH's horizon in the framework 
of a non-equilibrium statistics. Any further study, such as to choose a well-defined Hamiltonian, and to consider the 
possibility of creation or annihilation of geometrical particles, and so on, is very interesting. 

I close this paper with some comments. Due to the discreteness of spacetime itself at the Planck scale, there is 
a minimum length, namely the Planck length, which is similar to the lattice constant in condensed matter. Thus, 
mature methods and viewpoints developed from condensed matter physics and statistical physics can be used for 
reference in the study of the discrete quantum spacetime. In this direction, some groups have made interesting and 
enlightening efforts (see, for example, Refs. [9-12]). 
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APPENDIX A: PROOF OF EQ. (5) 

According to the definition (4) of (/ {Ak)) and the master equation (2), we have 

j t (/ (40> = £ ( ft 9 (A a ) J / {A).) £ £ g (A\) [-Wi (A, A\) P ({A} , t) 

{A} \a=l / i Xi 

+W t (Ai -» A 4 ) P ({A^ fc },i fc ,t)] 
= £(ns(40)/(40 £ £.9(i») [-W^A^A^PUA}^) 

{A} \a=l ) i(i^k) A t 
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+Wi (Ai^ A^p([A^ k },A k ,t) 

+ E ( II 9 (40 J / (40 E 5 (A) [-w* - 4) P ({^} - *) 

+W k (A k -» A fe ) P ({A j¥fe } , i fe , ) 



(Al) 



Looking at the i ^ fc term of (Al), 



N 



( ? ;^fc)tcrm= J] J] a(A a ) ] /(A fe ) 



x E 

i(»^fe) 



^ g {Ai) g (Ai) W % [A, - i 4 ) P ({A} , t) 
+ ^ g (A,) g (A^ Wi (Ai -> A<) P ({A,y fc } , i fe , i) 



(A2) 



it is easy to see that this term equals to zero, as long as Ai exchange with Ai before doing sum for Ai and Ai. So the 
surplus term of (Al) is only the last term (i = k): 



N 



Jt (/ (40) = E II 9 (A a ) J (A k ) £ g (i fe ) [-W k (M ^ A k ) P {{A} , t) 

{A} \a=l / i fc 

+Wfe (ife -» Afc) P ({A j¥ fe},ife,i) 



= -E (il^)) /(4) (^g(i fe )^(A fc ^i fe ) ] P({A},i) 
+ E <7 04i)---ff(A fc ) ff (i fc ) ••• 5 (A JV )/(A fc ) 
■xW k (A k -» A fe ) p({A j¥fe },i fc ,t) 

= - (/ (A fc )) + e ( n ^ (4*) ) ( e 5 (4) / (4) ^ (4 - a) ] p ({A> , *) 

in which, the normalized condition Y^A k 9 (^ fe ) ^ k (4 ~ * 4) = ^ an< ^ ^ e technique of exchange of ij for Ai were 
used. Hitherto, Eq.(5) have been proven exactly. 
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